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Motivated by a search for experimental probes to access the physics of fractionalized excitations 
called spinons in a spin liquids, we study the interaction of spinons with lattice vibrations. We 
consider the case of algebraic spin liquid, when spinons have fermionic statistics and a Dirac-like 
dispersion. We establish the general procedure for deriving spinon-phonon interactions which is 
based on a symmetry considerations. The procedure is illustrated for four different algebraic spin 
liquids: 7r-flux and staggered-flux phases on a square lattice, 7r-flux phase on a kagome lattice, and 
zero flux phase on a honeycomb lattice. Although the low energy description is similar for all these 
phases, different underlying symmetry group leads to a distinct form of spinon-phonon interaction 
Hamiltonian. The explicit form of the spinon-phonon interaction is used to estimate the attenuation 
of ultrasound in an algebraic spin liquid. The prospectives of the sound attenuation as probe of 
spinons are discussed. 
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I. INTRODUCTION 

Unambiguous experimental identification of a spin liq- 
uid, 1 an exotic ground state of a spin system in a di- 
mension larger than one without a magnetic order re- 
mains an open question. 2 A number of theoretical sce- 
narios leads to a ground state with charge-neutral exci- 
tations, which carry spin-1/2 quantum number and have 
fermionic statistics. These excitations, called spinons, 
may have a Fermi surface or a Dirac spectrum, 2,3 and 
are usually strongly coupled to a gauge field. The lat- 
ter case, so-called Dirac spin liquid phase, will be the 
primary subject of attention in the present paper. 

Naively, one would expect that the presence of 
fermionic excitations in the system could be easily tested 
experimentally. However, measurements of different 
thermodynamic quantities, such as spin susceptibility or 
specific heat, often require subtraction and extrapolation 
to zero temperature, which can be ambiguous. Transport 
measurements arc limited to a heat conductivity due to 
neutral character of spinons. The heat transport mea- 
surements are difficult to perform at low temperature. 
Finally, neutron scattering, potentially a direct probe of 
spinons, 4 ' 5 requires large single crystals which are not 
always available. The difficulties with experimental de- 
tection of spin liquid phases has lead to a number of 
theoretical proposals. These include but are not lim- 
ited to Raman scattering, 6 inelastic X-rays scattering, 7 
Friedel oscillations, 8 electron spin resonance, 9 impurity 
physics, 10-12 and optical conductivity. 13 ' 14 

Coupling between spinons and phonons can open an- 
other channel of decay for phonons, and change the at- 
tenuation of ultrasound. Thus, the sound attenuation is 
another potential probe of spinons. 15 The case of spin 
liquid with spinon Fermi surface has been recently con- 
sidered by Zhou and one of us in Ref. 16. The spinon- 
phonon interaction in the long wavelength limit was de- 
duced from hydrodynamical arguments. 16-18 Assuming 
that electrons stay in equilibrium with lattice (due to 



presence of impurities) and making canonical transfor- 
mation to the moving frame one can easily derive inter- 
action Hamiltonian. The same interaction Hamiltonian 
can be reproduced from microscopical considerations us- 
ing so-called deformable ions model. 19,20 

However, as was pointed out in Ref. 16, spinons with 
a Dirac spectrum require a different treatment. Here 
we address this problem and consider the contribution 
of spinons to the attenuation of ultrasound in a Dirac 
spin liquid. In contrary to the case of electrons 18-20 or 
spinons with Fermi surface, 16 it appears that there is no 
universal form for the interaction of spinons with acous- 
tic phonons in Dirac spin liquid. This is related to the 
spinor nature of spinons with Dirac dispersion. Simi- 
larly to graphene, spinor structure describes the charac- 
ter of wave function on different sublattices, i.e. on the 
microscopic scale. Consequently, one possible route of 
deducing the form of spinon-phonon interaction would 
be to find the change of the Hamiltonian of electronic 
subsystem induced by the long wave modulation of the 
lattice parameters, starting from a microscopic Hamilto- 
nian. This procedure, although giving explicit values of 
coupling constants is not universal. It depends on the mi- 
croscopic implementation of a different spin liquid phase. 
Moreover, it is difficult to guarantee that one finds all 
possible terms in the interaction Hamiltonian. 

We adopt a different approach and use symmetry con- 
siderations to find a spinon-phonon interaction Hamil- 
tonian. Similar route has been recently used to deduce 
electron-phonon interaction in graphene. 21 ' 22 There, rep- 
resentations of the lattice symmetry group of the hon- 
eycomb lattice on continuous Dirac fields were used to 
find all possible symmetry allowed electron-phonon cou- 
plings. 22 However, in a Dirac spin liquid, also referred 
to as an algebraic spin liquid, the notion of symmetry 
group has to be extended to the projective symmetry 
group. 23 ' 24 

In this paper we generalize the derivation of spinon- 
phonon interaction Hamiltonian to the case of projective 
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representation of lattice symmetry group. The proce- 
dure is straightforward, and it requires studying the rep- 
resentation of symmetry group on spin-singlet fermionic 
bilinears. We consider four different realizations of the 
Dirac spin liquid: 7r-flux 25 and staggered-fiux 26 phases 
on a square lattice, as well as 7r-flux phase on a kagome 
lattice 27 and a Dirac spin liquid phase on a honeycomb 
lattice. Within low energy effective field theory, all these 
phases can be described in terms of a Dirac excitations, 
coupled to a gauge field. 25-31 Nevertheless, these phases 
retain the information about their microscopic origin. 
This information is encoded in the projective represen- 
tation of spinon operators under the action of the corre- 
sponding symmetry group. 

As we will see below, the projective character of the 
representation of symmetry group has a profound con- 
sequences on a spinon-phonon interaction. For all alge- 
braic spin liquid phases considered here except for the 
Dirac spin liquid on a honeycomb lattice, only the cou- 
pling to the density of spinons is allowed by symmetry 
at the leading order. However, we find that this coupling 
is screened due to the gauge field. The sound attenua- 
tion coming from the next order terms is suppressed by 
an extra factor and behaves as T 3 at low temperatures. 
In contrast, the Dirac spin liquid on a honeycomb lattice 
has a sound attenuation oc T. Although the sound atten- 
uation in all cases is suppressed compared to the case of a 
spin liquid with a Fermi surface, the spinon contribution 
is still the dominant process at low temperature and ex- 
perimental observation of the attenuation of ultrasound 
due to spinons may be possible. 

The paper is organized as follows. Section II intro- 
duces the low energy description of Dirac spin liquid and 
projective symmetry group. We use the 7r-flux phase on 
a square lattice as an example. Next, we discuss interac- 
tion of spinons with phonons in Section III. We describe 
the general procedure of obtaining spinon-phonon inter- 
action Hamiltonian from symmetry considerations. It is 
later illustrated in more details for the 7r-flux phase on 
a square lattice. Having found the interaction Hamilto- 
nian, in Sec. IV we study the sound attenuation, con- 
centrating on attenuation of longitudinal sound. Finally, 
in Sec. V we summarize and discuss our results. Basic 
facts from the representation theory of finite groups and 
details on the calculation of sound attenuation are given 
in Appendix. 



II. LOW ENERGY DESCRIPTION OF 
ALGEBRAIC SPIN LIQUIDS 

We review the description of algebraic spin liquid 
fixed point using language of low energy effective field 
theory. 25-31 This description suits our purposes since we 
are interested in coupling between acoustic phonons in 
the low energy limit. In addition, it provides universal 
framework applicable to a variety of different algebraic 
spin liquid phases. In short, three ingredients are needed 



in order to specify low energy field description of a given 
algebraic spin liquid phase. These are low energy Hamil- 
tonian for continuous fields, gauge group and representa- 
tion of projective symmetry group specified by its action 
on fields. Below we review these ingredients for the gen- 
eral case, as well as illustrate them for the 7r-flux phase 
on a square lattice (irFn phase). We do not list details 
for the staggered flux phase on a square lattice (sFn), 
7r-flux on kagome lattice (7rFo), and Dirac spin liquid 
phase on honeycomb lattice (OFo). The reader is re- 
ferred to Refs. 25-27 for more details, specific for these 
phases. 

A. Effective field theory, gauge group and 
projective symmetry group 

The starting point is the spin 5 = 1/2 model on some 
(not necessary Bravais) two-dimensional lattice. The 
symmetries of spin Hamiltonian are assumed to include 
SU(2) spin rotations, time reversal and the full lattice 
group of a given lattice. We write Hamiltonian as 

H = J2JijSi-S j + ..., (1) 

(ij) 

where bold indices i,j denote lattice sites and sum goes 
over nearest neighbor pairs of sites. Ellipses denotes 
other short-range interaction terms required to stabilize 
required spin liquid phase. 

In order to get access to phases with no spin order, we 
use slave-fermion mean field theory. We represent spin 
operators using spinon operators /i jQ , a =t, 4- : 

S i = ^fLvapfjp- (2) 

The mapping between Hilbert spaces is exact, pro- 
vided one imposes a constraint of no double occupancy, 
fiafia = 1 7 where summation over repeating indices is 
implied. Such representation of spin has SU(2) gauge 
redundancy. Therefore, mean field decoupling of spin in- 
teraction has to include SU(2) gauge field on the links 
of the lattice. Saddle points of the mean field theory, 
depending on their structure 23 ' 24 may break this SU(2) 
symmetry down to U(l) or Z2. 

We will be mostly interested in phases with U(l) gauge 
symmetry (the only SU(2)-symmetric algebraic spin liq- 
uid phase considered here is the ttFo phase 24 ). General 
lattice gauge theory Hamiltonian at the saddle point is 

(ij) plaq. 

+ E^ e ^ i /i Q / J o+H. c], (3) 

(ij) 

where we assume that the choice of Xij leads to a Dirac 
spectrum. The gauge redundancy lead to appearance of 



dij, a compact U(l) vector potential living on the bonds 
of the lattice, and eij which is canonically conjugate elec- 
tric field taking integer values. We also defined </>□, liv- 
ing on the dual lattice, as the lattice version of curl of 
gauge field a^-. Gauge constraint is (div e)i+fj a fia = 1- 
Although initially vector potential does not have any ki- 
netic terms, these will be generated due to coupling with 
fcrmions and are written in (3) with coefficient K. For 
K = 0, h^> \xij\> one recovers spin model from Eq. (3). 
Mean field is a good approximation when K ^> \xij | ^> h 
and describes the algebraic spin liquid phase. It has been 
argued that this phase is stable. 25,26 Monopoles are irrel- 
evant, leading to the gauge group becoming non-compact 
in the low energy limit. Therefore it should be accessible 
starting from the spin Hamiltonian in Eq. (f) for some 
range of initial parameters. More detailed discussion of 
stability of algebraic spin liquid phase is presented in 
Refs. 25 and 26. 

Provided that choice of Xij leads to a Dirac spectrum, 
we use continuous fermionic fields to write low-energy 
Hamiltonian. In the momentum space it reads 



H = V F 



d 2 k 



(4) 



The indices a, a and a in continuous eight-component 
fermion field ipaaa label spin, Dirac valley and sublat- 
tice respectively. In what follows, we will use three dif- 
ferent sets of Pauli matrices acting in different spaces. 
Pauli matrices acting on the spin indices are denoted as 
{a x , a y ,<T z }. The second index distinguishes between dif- 
ferent Dirac points (valleys) and we label Pauli matrices 
acting in this space as {[i x , fi v , /i z }. Finally, Eq. (4) in- 
volves Pauli matrices {t x ,t v ,t z } acting on spinor (sub- 
lattice) indices. Below we will omit the sign of the tensor 
product implying e.g. a v r x — a v <E) 1 <8> t x . 

A particular choice of anzatz Xij naively violates some 
symmetries of the original Hamiltonian. However, repre- 
sentation Eq. (2) is invariant under the gauge transfor- 
mations 



fi he 



iOi 



(5) 



Consequently, there is a freedom in the choice of the ac- 
tion of different symmetries: one can supplement them 
with some gauge transformation. Requiring mean field 
Hamiltonian (3) to remain invariant we can fix this free- 
dom and show that all symmetries of original spin Hamil- 
tonian remain unbroken. 

The fact that the action of symmetries on fermionic 
fields is supplemented by gauge transform has a deep 
consequences. Rather than forming usual representation 
of lattice symmetry group, fermions are said to realize 
•projective representation of lattice symmetry group. 23 ' 2 
This representation is fully specified by the action of 
the lattice symmetry group generators on the fermionic 
fields. Knowing the mapping from the lattice fields fi a 
to ip, one can easily find the action of generators on the 
continuous fermionic field. Below, we will demonstrate 
this procedure for the 7rFn phase. 




FIG. 1. Choice of the anzats for the ttFo phase. Red dashed 
line encloses the unit cell. Numbers indicate the labeling of 
sites, used in the main text. Hopping in the direction of 
arrows is proportional to +i. 



B. Example: 7r-flux spin liquid on a square lattice 

We use the the 7rFn phase 25 to illustrate the abstract 
construction presented above. The starting point is the 
choice specific values of Xij m Eq- (3) to fix the mean field 
anzats. Our notations are different from those in Ref. 25. 
We quadruple our unit cell, so it includes four sites. Sites 
within the new unit cell are labeled by an index m = 
1 ... 4 as in Fig. 1. The Xij an d a° are chosen as 



Xi 



= i, 



xu+y = (-)M; og = o. (6) 



Resulting mean field Hamiltonian with the spin index 
and gauge field omitted is: 

H °* = "^X>(/rt + fldfrl + KfU + ft+a 2 l)fr4 
r 

K/r+aa _ fr3)fr2 + Kfr+aiA ~ fr+a 2 2)fr3 + H. C.] , 

(7) 

where Bravais lattice vector r = n\ai +rt2<i2 with m, n2 
integers labels unit cells, and = 2x, a 2 = 2y. Making 
Fourier transform 



fmi a/JV ^ ^ ^ fkmi 



(8) 



where TV is the number of unit cells, we obtain the Hamil- 
tonian in momentum space 



HmS — o fkm-Hmn{k)fkn, 



(9) 



where 4x4 matrix H(k) is given by 



/ -1 + K{ 



H mn (k) — i 



1-Kt 


\l-K 2 









1 -K* 



-l + K* 





l + K 2 1-Ki 
-1 + Kf 



(10) 

with the notation K\ t2 — e 1 1,a . The momentum is 
measured in the units of the inverse lattice constant of the 
original lattice (i.e. before enlarging the unit cell), which 
is set to one in what follows. We choose the Brillouin 
zone as k x ,k y € [— tt/2, 7t/2). The energy levels of the 
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Hamiltonian (9)-(10) are doubly degenerate. Therefore 
there is only one doubly degenerate gapless Fermi point 
in the Brillouin zone, that is located at the momentum 
Q = (0, 0). In the vicinity of this point the dispersion is 
described by Dirac fermions and we can define continuum 
fermion fields. We choose the following representation for 
two copies of spinor field ipaa{k): 



1 ( l/fc2 + fk4 
y/2 \— i/fel — /fe3 J ' 

1 ( fki + i/fc3 
y/2 \~fk2 — i/fe4 



(11a) 
(lib) 



where the index a = 1,2 labels different Dirac points, 
K±, both located at the T point of the Brillouin zone. 
Consequently, the low-energy continuum Hamiltonian is 
written as 



cPk 

{2-Kf 



^aa(k X T X + k y T y )i> a 



(12) 



where made spinor index implicit. Fermi velocity, pro- 
vided lattice constant is set to one, coincides with J, 
Vp = J. 

The continuous fields t/> aa (k) realize a projective repre- 
sentation of the lattice symmetry group. This represen- 
tation is fully specified if one knows the action of group 
generators on continuous fields. The relevant lattice sym- 
metry group in this case, denoted as C' 4v , is different from 
the point symmetry group of square C& v due to enlarged 
unit cell. Consequently, the group C^ v , in addition to 
elements from C^ v , contains translations by a unit vec- 
tors of square lattice along x and j/-axes. The full set of 
generators, action of which is to be specified below, con- 
tains rotation for angle n/2, TZ^/2, reflection of x-axis, 
TZ X (these are generators of Ci v ) and lattice translation 
by vector Oi, T x . In addition we also have to specify the 
action of time-reversal symmetry, T, and charge conju- 
gation operators, C. 

Let us illustrate the derivation of the action of reflec- 
tion generator on the continuous fields. Action of 1Z X on 
the unit cell may be symbolically shown as 



n x 



f ; 




< 
< 





(13) 



From here we can understand the action of 1Z X on the 
lattice fermion fields, as exchanging fermionic operators 
with indices 1 f-> 2 and 3 f-> 4. However, in order to 
leave Hamiltonian (7) invariant, this has to be supple- 
mented by gauge transform that changes the sign of all 
hoppings [thus reversing the direction of all arrows in the 
r.h.s. of Eq. (13)]. One can easily check that transforma- 
tion 

frl fr'2, fr2 ~» ~fr'l, (14a) 
/r-3 fr'4, frA -> -fr'3, (14b) 

r' = K x r = (-r x ,r y ) (14c) 



leaves Hamiltonian (7) invariant. The action of TZ X on 
operators f rn is easy to translate into representation of 
1Z X on the continuous fermionic fields using their defi- 
nition (11). In terms of different sets of Pauli matrices 
introduced above, it can be written as 



ip K x ip with K x = \[i z T y 



(15) 



Using analogous procedure we get the representation of 
the remaining generators to be 



T x = [i v . 



(16) 
(17) 



Finally, for the 7r-flux phase there exist two additional 
SU(2) transformations not related to lattice symmetry 
group. These are time reversal, T, defined as a antiu- 
nitary operator which flips the direction of spin opera- 
tor (2), and charge conjugation, C. The latter operation 
may be viewed an 577(2) rotation in the spin space by it 
around y-axis, 25 supplemented by the particle-hole trans- 
formation. On the lattice level (in momentum space) 
these can be conveniently represented by 



T : fkal,3 «-» /feo-1,3) fka2A ^ ~ f ko2 



C ■ fk^n I f—htn' fk^n ^ /. 



■fefn' 



(18a) 
(18b) 



where T-reversal also includes complex conjugation op- 
eration, spin indices were restored. Although these defi- 
nitions may look counterintuitive, one can check that the 
charge conjugation (18b) indeed leaves the spin operator 
invariant, whereas the time-reversal symmetry, defined 
as in (18a), flips the direction of spin operator in Eq. (2). 
Mapping this action to continuous fermionic fields, we 
have: 



T: 
C : 



(19a) 
(19b) 



As we pointed out earlier, the representation of the lat- 
tice symmetry group on fermions is projective. This can 
be easily seen from the action of generators, Eqs. (15)- 
(17), if one tries to test some group identities. For exam- 
ple, (lZn/2) 4 is a trivial transformation. However, using 
the explicit form of the representation of H n /2 for con- 
tinuous fermionic fields, Eq. (16), we find 



(ft* 



/2j 



1. 



(20) 



Thus all group identities hold only up to some gauge 
transformation, which leaves the Hamiltonian invariant. 



III. SPINON-PHONON INTERACTION 

Having a low energy description of Dirac spin liquid 
phases at our disposal, in this section we consider the 
spinon-phonon interaction. 

As we explained in the introduction, the hydrodynamic 
approach, 16-18 applicable for the case of spin liquid with 
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the Fermi surface, is not straightforward to use in our 
case. It is the presence of spinor structure, inherently re- 
lated to the microscopic details such as two inequivalent 
sublattices, that prevents application of hydrodynamical 
arguments. Of course one can always resort to the mi- 
croscopic derivation of spinon-phonon interaction. De- 
spite the advantage of giving specific values of coupling 
constants, this route is highly non-universal and is not 
guaranteed to yield all possible couplings. 

We present universal procedure for finding all possible 
terms in spinon-phonon interaction Hamiltonian, allowed 
by symmetry. First, we introduce phonons and the gen- 
eral form of the spinon-phonon interaction Hamiltonian. 
After this the general idea behind the procedure is de- 
scribed. Implementation of this procedure for the 7rFn 
phase serves as an example. Finally, we present results 
for other Dirac spin liquid phases and discuss the under- 
lying physics. The derivation of these results extensively 
relies on a representation theory for finite groups. Neces- 
sary concepts, as well as basic facts about point groups 
of square, kagome and honeycomb lattices are listed in 
Appendix A. 

A. Spinon-phonon interaction Hamiltonian from 
symmetry considerations 

We start with specifying conventions for the spinon- 
phonon interaction Hamiltonian. It is written using the 
operator jy s _ p h(fc, q) as 

H s . ph = ^( fe + q)H s . ph (k, g)^(fc). (21) 

In what follows we will refer to the operator i? s _ p h(fc, q) 
itself as a spinon-phonon interaction Hamiltonian. Nor- 
mally, the operator iT s _ p h(fc, q) obtained from the pro- 
cedure described above, would contain only zeroth or- 
der terms in the distance from the Dirac point, fc, 
H s - P h{k, q) — -ffs-ph(o)(?)- As we shall see, in some cases, 
all such terms vanish. Then, to find a non-zero interac- 
tion Hamiltonian, we allow the presence of terms, linear 
in k, -ff s _ p h(i) (k, q) and the total Hamiltonian will be 
written as a sum: 

H s - ph (k,q) = i? s _ ph (o)(q) + H s _ ph(1) (k,q), (22) 

In what follows, we restrict ourselves to the first 
non-vanishing term in this expansion. Phonons enter 
if s _ p h(fc, q) via the q Fourier component of displacement 
field, u(t,r). In the second quantized language the dis- 
placement field is written as 

(23) 

where index \x — L,T labels longitudinal and transverse 
modes of acoustic phonons, and e qfi is the correspond- 
ing polarization vector. The dispersion of phonons is 



assumed in the form u! q — v s \q\, where v s is the sound 
velocity. The p is defined as a mass density per layer, and 
S is the area. For simplicity we consider only in-plane 
phonon modes. 

Although we work in continuum limit, it is the lat- 
tice symmetry group and its representations, which de- 
termines the properties of low energy (acoustic) phonons 
and spinon excitations. Phonons are described using vec- 
tor u(t, r), describing displacement at a given point r due 
to deformation. As a uniform displacement of the entire 
lattice, u(t,r) = u , leaves system invariant, acoustic 
phonons can couple to spinons only via spatial deriva- 
tives of u(r) [we ignore coupling to the time derivative, 
as it is suppressed by the ratio of sound and Fermi veloc- 
ities]. Set of all spatial derivatives, diUj(r) or —iqiUj(q) 
in the Fourier space, transforms as a rank-two tensor un- 
der lattice symmetry group. Representation of a lattice 
symmetry group on a rank-two tensor can be split as a 
sum of irreducible representations. Symbolically this is 
written as 

EixE 1= J2 ® D f , ( 24 ) 

3 

where E\ is vector representation, and D? h are (possibly 
repeating) irreducible representations. Acoustic phonon 
modes can be classified using irreducible components, 
present in this decomposition. 

Spinons have fermionic statistic, thus minimal cou- 
pling to phonons must involve bilinears of ip field. In 
contrary to phonons, continuous spinon fields ij) realize 
projective representation of the lattice symmetry group. 
Action of lattice symmetries on i/j in general includes the 
gauge transformation, and all identities between genera- 
tors are valid modulus gauge transformation [for exam- 
ple, see Eq. (20)]. Similar to a single field tp, general 
spinon bilinear also realizes projective representation of 
lattice symmetry group. However, there exists a subset 
of spinon bilinears which transform under regular rep- 
resentation of symmetry group. For the case when the 
gauge group is SU(2) these are bilinears which are sin- 
glets under SU(2). Whereas for abelian gauge groups, 
like U(l) or Z2, all bilinears realize regular representa- 
tion, as gauge component cancels. 

Regular representation of the lattice symmetry group 
on (a subset of) spinon bilinears can be split into irre- 
ducible representations Djf , 

G^=J2®Df*- (25) 
3 

We note, that invariant fermionic bilinears for 7rFn, sFn 
and 7rF§s phases were identified in Refs. 25-27. This 
corresponds to finding all trivial components contained 
within decomposition (25). 

The product of two irreducible representations, Df h x 
Dj' contains a trivial representation within itself if, 
and only if these representations coincide, Df h = ^ . 
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Representation Ax A 2 B[ B' 2 Ex E 2 E{ E s Ex E 5 

Basis 1 t z li z n z r z t x ,t v fi x , fi y ^r x ,^r y ^ x r z ,Li y r z ^ T y , fJ, y T x fj, x r x ,Li y T y 

T-inv ---- + + + + 

C-inv - + +--- + + 



TABLE I. Explicit form of basis in terms of tensor products of Pauli matrices for irreducible representations of C' iv contained 
within G^f Last two rows show properties of basis elements under time-reversal and charge conjugation. Plus implies 
invariance, whereas minus indicates a change of sign under the action of corresponding symmetry. 



As the spinon-phonon interaction Hamiltonian has to be 
invariant under the action of the symmetry group, we 
can construct all symmetry allowed couplings by pairing 
identical irreducible components between splittings (24) 
and (25). The presence of additional symmetry opera- 
tions, such as time-reversal or charge conjugation may 
impose further restrictions on the obtained set. 

B. Example: derivation for the 7r-flux phase 

We use the irF phase as an example for an illus- 
tration of the abstract procedure outlined above. For 
phonons, the underlying symmetry group is C± v (see Ap- 
pendix A 2). Using characters Table V, we find explicit 
form of the decomposition (24) for the present case: 

Ei x Ex = A 1 ®A 2 ®B 1 ®B 2 . (26) 

Ax here and in what follows always denotes the trivial 
representation. All other representations are also one- 
dimensional, thus the action of corresponding group ele- 
ments can be inferred from Table V in the Appendix. In 
terms of components of two vectors (q x , q y ) and (u x ,u y ), 
the basis functions of these representations are 

A\ '■ u xx + Uy y) A 2 : u xy — u yx , (27a) 
Bx-u xx -u yy , B 2 :u xy + u yx , (27b) 

where we introduced shorthand notation 

Uij = QiUj. (28) 

While introducing the ansatz for 7rFn phase in Sec- 
tion II B we used unit cell consisting of four lattice sites. 
This allowed us to write relations between continuous 
fields and microscopic spinon operators in a simple form. 
However, the price to pay is that spinons now trans- 
form under the symmetry group C' iv which is larger 
than point symmetry group of the square lattice. In ad- 
dition to transformations from the point group of the 
square Ci V , group C' iv includes lattice translations by 
unit vector in x and y directions. The details about ir- 
reducible representations of group C' iv are worked out in 
Appendix A 2. It has eight one-dimensional irreducible 
representations labelled as Ax,2i Bx, 2 and A\ 2 , B[ 2 in Ta- 
ble VII, and six two-dimensional representations denoted 
as Ex, E[, £^2,...,5- To find the splitting of representation 



of C' iv on spinon bilinears into irreducible components we 
use natural basis: all spin singlet bilinears can be enu- 
merated using tensor products of Pauli matrices acting 
in sublattice and valley space, 

V^MV-, Me{l,r ! ,/i ! ,rV}. (29) 

With the help of the characters Table VII, the 16- 
dimensional representation G^f^ is reduced into direct 
sum of four one-dimensional and six two-dimensional rep- 
resentations as 

G>t^ = M © A 2 © B'x © B' 2 © Ex © E[ ® E 2 

®E 3 ®E 4 ®E 5 . (30) 

Using the explicit basis (29), we can find to what irre- 
ducible representation a given matrix belongs. Identity 
matrix 1 corresponds to the trivial representation, A\. 
Next, one can check that both matrices n z , t z , and their 
product, [i z t z , are invariant (up to a sign) under the ac- 
tion of all generators of C' 4v , Eqs. (15)-(17). Therefore 
the matrices /i z , t 2 , and [i z t z form the basis of one- 
dimensional representations A 2 , B[ and B' 2 respectively. 
Matrices (T x ,T y ), ([i z t x , [i z t v ), and (fj, x ,fx v ) constitute 
basis of two-dimensional irreducible representations Ex, 
E[ and E 2 respectively. Finally, after some algebra, the 
remaining six matrices from (29) can be split into pairs 
that realize the basis for representations Ea,...,E& as 
shown in the Table I. The last two rows in Table I display 
the symmetry of corresponding matrices under the action 
of time reversal and charge conjugation operations. 

Explicit decompositions, Eqs. (26) and (30) give us 
allowed couplings between spinons and phonons. Only 
identical irreducible representations can be coupled be- 
tween themselves. Comparing Eqs. (26) and (30) we see 
that only two first terms in both direct sums coincide. 
Thus one may expect the allowed couplings to be de- 
scribed by contraction between Ax {A 2 ) components in 
different sums. However, according to Table I, the bilin- 
ear tp'^ip is odd under both time-reversal and charge con- 
jugation, and if, ! >T z ip is odd under time reversal. As differ- 
ent components of invariant under time-reversal 
and charge conjugation, we conclude that at the leading 
order no couplings of spinons to phonons are allowed by 
symmetry, Hl^ {0) (q) = 0. 

To find a non-zero coupling of spinons to phonons, we 
allow for the presence of spinon momentum, fe in the cou- 
pling Hamiltonian. This corresponds to the next order 



7 



Representation A 2 Bi B[ A' 2 E[ E 3 E 4 E 5 

Basis t z r z a z fi* 1 t*' v ,H*t x » , r z ^ x ' y ^^(t^t") (i x »(t x ±t v ) 

T-inv - - - - + + 



TABLE II. Explicit form of basis for different irreducible representations of C± v contained within G B J^. Action of the group 
generators coincides with Ref. 26. Last row summarizes the transformation of basis elements under time-reversal symmetry. 



in the expansion around the Dirac points. The spinon 
momentum, fc, transforms under the usual vector rep- 
resentation Ei [this can be inferred from the fact that 
Dirac Hamiltonian, Eq. (12) is invariant] and is invariant 
under time reversal, but odd under charge conjugation. 
The product Ei x G^^ is reduced as 

E l x Gf H , = A l ®A 2 ®B l ®B 2 ® A[ ® A' 2 ® B[ ® B' 2 
®2(E 1 ®E' 1 ®E 2 ®E 3 ®E i ®E s ). (31) 

The first four irreducible representations, which are of 
interest for us [cf. with Eq. (26)], originate from the E\ 
component within G^^. Consequently, their basis is 
analogous to Eq. (27). The first four irreducible represen- 
tations A\ . . . B 2 from Eq. (31) can be coupled to corre- 
sponding irreducible representations in Eq. (26). For ex- 
ample, coupling representation A\ with basis k x r x + k y T v 
to A\ component with basis u xx + u yy results in contri- 
bution 

= (««X + Uyy) [k x T* + ky T y] , (32) 

with a phenomenological coupling constant g j^ . Depen- 
dence on fc, q will be suppressed for brevity in what fol- 
lows, i? s _ph(i) = H s . p ui-)(k, q). Collecting all contribu- 
tions at this order, and rearranging phenomenological 
coupling constants (e.g. g^\ = g^±g^) we get the most 
general form of the spinon-phonon interaction Hamilto- 
nian in the 7rF phase to be 

-^S^ph = Si \ U Xxk X T X + UyykyT V ) 

^ 9 2 \ u xyk x T + Uy X kyT^) + g 3 (uy X k x T + u X ykyT y ) 

94 {u xx kyT^ -\~ Uyyk x T 

')■ (33) 

C. Results for the sFn, 7rFs£s, and OFo phases 

After detailed derivation of spinon-phonon interaction 
for the 7r-flux phase on a square lattice, we present results 
for other algebraic spin liquid phases considered in this 
work. 

The staggered- flux phase on a square lattice is simi- 
lar to the 7rFn phase, considered above. However, in 
contrary to the 7r-flux phase, there is no charge conjuga- 
tion present among additional symmetries. As we shall 
shortly demonstrate, due to reduced symmetry, the num- 
ber of allowed couplings is going to be larger. Splitting 



of G^f° into irreducible components works as 

Gffij, = A 2 ®B 1 ®B' l ®A! 2 ® 2E[ ® 2E 3 ®E 4 ® E 5 . 

(34) 

The basis of corresponding components in terms of tensor 
product of Pauli matrices are listed in Table II. One can 
notice, that all one-dimensional representation present 
in the decomposition (34) are odd under time-reversal. 
Just like the case of the 7rFn phase, no couplings with 
phonons are allowed by symmetry at this order. To find 
non-zero coupling we consider next order in fc. 

Explicit expression for -ff s _ p h(i) (fc) can be found us- 
ing decomposition of the product E[ x [where E[ 
corresponds to spinon momentum 32 ] into irreducible rep- 
resentations, 

E 'i x Gf,° = 2(Ai ®A 2 ®B 1 ®B 2 

®E 1 ®E[®E 2 ®E 3 ®E i ®E 5 ). (35) 

The overall factor of two indicates that there are two 
distinct copies of each irreducible representation in the 
decomposition. Four one dimensional irreducible repre- 
sentation in the first line of Eq. (35) coincide with the de- 
composition of Ei x Ei into irreducible representations. 
As all irreducible components are encountered twice in 
the decomposition (35), we have eight different couplings 
between spinons and phonons. 

The basis for the two copies one-dimensional represen- 
tations Ai . . . B 2 in Eq. (35) can be deduced using the 
fact that all these irreducible representations originate 
from the tensor product E[ x E[ : 

A x : k x r x + k y T y , A 2 : k x r y - k y r x (36a) 
Bi : k x r x - k y T v , B 2 : k x r y + k y r x . (36b) 

The bases for the second copy of irreducible represen- 
tations in Eq. (35) have the same form as in Eq. (36), 



Representation A\ 


A 2 Ei Fi 


F 2 


F3&F4 


Basis 1 


r z t x ' v r z a x ' v ' z 




x,v , x,y,z 


T-inv + 


- - + 




+ 



TABLE III. Explicit form of basis for irreducible representa- 
tions of C' 6v contained within G^J®. Notations and action of 

group generators coincide with those used by Hermele et. al?' 



Representation A\ 


Bi A2 B2 Ei 


E 2 


Ai 


Bi 


A2 B2 Ei 


E 2 


Basis 1 


H z r z fj, z T z r x ' y 




^r z 


fI V T Z 


IX* n y n x r y ,-fj, x r x 


H v t x ' v 


T-inv + 


- - + - 


+ 


+ 


+ 


- - + 


+ 



TABLE IV. Irreducible representations of Ca v contained within G^ t ^ and their basis. Each irreducible component occurs twice: 
first six representations in the Table are diagonal in the valley space, whereas remaining six are their off-diagonal counterparts. 
Adopted from Table III in Ref. 22. 



but with an extra Pauli matrix fj, z . Physically, this cor- 
responds to the fact that the anisotropy in Fermi veloc- 
ity is not prohibited by symmetry in the sFd phase. 26 
Time-reversal invariance does not reduce the number of 
allowed couplings, as all combinations in Eq. (36) are 
invariant under T . The resulting spinon-phonon interac- 
tion Hamiltonian, has eight terms, four of which coincide 
with the case of the 7rFn phase, Eq. (33), while remain- 
ing four terms contain an extra \x z and correspond to a 
phonon-induced valley anisotropy in a Fermi velocity. 

The n-flux phase on a kagome lattice has a symmetry 
group C' ev . Just like in the case of square lattice, it is an 
extension of a conventional symmetry group of hexagon, 
Cq v , to the case of enlarged unit cell. 27 The point group 
of hexagon governs the properties of phonons. The gener- 
ators of Cq v are rotations for an angle of 7r/3, and 
reflection with respect to cc-axis, 1Z X . Using character Ta- 
ble VIII, the tensor product of two vector representations 
is decomposed as 

E x x E x = A x ® A 2 ® E 2 . (37) 

On the other hand, using character table for the group 
C' 6v from Ref. 27, one can reduce representation on 
fermion bilinears as 

g$* = A x ® a 2 ® Ex ® f x ® f 2 ® f 3 ® f 4 . (38) 

Comparing Eqs. (37) and (38), we see that, in principle, 
couplings between corresponding A x and A 2 components 
are possible. However, if we consider properties of rep- 
resentations under time reversal, listed in Table III, we 
see that basis of A 2 is odd under time reversal and only 
coupling to density fluctuations remains: 

HZlt(0)=9^(u xx + u yy )l. (39) 

This coupling vanishes for transverse phonon modes, and 
we consider terms which are next order in k. Higher order 
terms can be found from decomposing of E x x G^J* into 
irreducible components, 

E i x = A i © ^2 © 2£d © E 2 © 2(F X © F 2 © F 3 © F 4 ), 

(40) 

and pairing those with irreducible representations con- 
tained within Ei x E\, Eq. (37). Only irreducible rep- 
resentations Ai, A 2 and E 2 in Eq. (40) are of interest 
as the same components are also present in decompo- 
sition (37). The basis for these representations can be 



written in terms of basis of E±, (k x ,k y ) and E\ compo- 
nent within G*%*, [t = ( T x ,T y ), see Table III]: 

A x : k ■ t (41a) 
A 2 : k x t (41b) 
E 2 : (k x r v + k y r x , k x r x - k v r y ). (41c) 

From here, we can read off the most general form of the 
spinon-phonon interaction Hamiltonian to be 

= 9 { i ] i u xx + u yy ){k ■ t) + g { 2 1] {u xy ~ u yx )(k x r) 

~\~93 ^ \_{ u xy+ u yx)(k x T V +kyT X )-\-(u xx —Uyy)(k x T X —k y T V )~^ . 

(42) 

The uniform phase on a honeycomb lattice is also gov- 
erned by the symmetry group C$ v . The notable differ- 
ence compared to the cases considered above, is that the 
honeycomb lattice is not Bravais and has a unit cell con- 
sisting of two atoms. Derivation of spinon-phonon inter- 
action is analogous to the case of graphene. 22 We neglect 
the optical phonon modes related to the presence of two 
atoms in the unit cell. Tensor product of two vector 
representations is given by Eq. (37). Whereas, decom- 
position of G°^ into irreducible representations works 
as: 22 

Gf f °=2(A 1 ®A 2 ®B 1 ®B 2 ®E 1 ®E 2 ), (43) 

where two copies describe matrices diagonal and non- 
diagonal in valley space. Bases of different irreducible 
components are given in Table IV. Terms which are off- 
diagonal in the valley space are not considered, as we do 
not allow for the intervalley scattering due to phonons. 
At the zeroth order in fc, in addition to the density cou- 
pling, Eq. (39), there are terms which do not vanish for 
transverse phonon modes, 

H S F ph = 91 [(UXX - Uyy)T X - {U X y + Uy X ) T V ] fl' . (44) 

Thus there is no need to consider next order in k. Note, 
that as basis of A 2 is odd under time reversal, we omit 
otherwise possible term (d x u y — d y u x )T z from Eq. (44). 

D. Comparison between different phases 

It is instructive to compare the above results for 
the spinon-phonon interaction in different realizations of 



Dirac spin liquid. In all derivations we considered Dirac 
fermions describing low energy excitations. The Dirac 
dispersion arises as an approximation of the band struc- 
ture in vicinity of K± points in the Brillouin zone. Con- 
sequently, the interaction with acoustic phonons may be 
understood from the influence of lattice deformations on 
the low energy band structure. The coupling of phonons 
to the density of spinons is very easy to explain from this 
perspective. The local changes in the volume of the lat- 



tice, described exactly by div u 



correspond 



to the density modulations of spinons, yielding the inter- 
action Hamiltonian (39). In the case of the OFo phase, 
remaining terms given by Eq. (44) are can be interpreted 
as a relative shift of K± points with respect to each other 
by lattice deformations. In other words, strain is trans- 
lated into a gauge field, which coupled with opposite sign 
in different valleys - well known effect for the case of 
graphene. 22 ' 33 

The presence of fluxes and non-trivial action of projec- 
tive symmetry group prohibits density coupling for ttFo 
and sFn phases. In the ttF& phase, the density coupling 
is the only allowed coupling at this order. To find a non- 
trivial couplings, we considered next order expansion in 
vicinity of the Dirac points. These couplings may be 
readily understood as a deformation of the band struc- 
ture in vicinity of K± point, which, nevertheless leaves 
the position of the Dirac point within the Brillouin zone 
intact. This is exactly what we see in couplings (33) and 
(42), which can be interpreted as the change in Fermi 
velocity, vp. Note, that the position of Dirac points in 
the Brillouin zone is non-universal, and depends on the 
choice of the implementation of the given phase. There- 
fore, it is natural, that the (physical) lattice deformation 
does not have any impact on the (unphysical) position of 
the Dirac points. 



IV. SOUND ATTENUATION 

We continue with a discussion of observable conse- 
quences of spinon-phonon interaction. Interaction of 
acoustic phonons with gapless spinons opens another 
channel for decay of phonons. Thus, it is expected to 
contribute to the attenuation of ultrasound. To get an es- 
timate of the this effect, we perform a simple calculation 
in this Section. As an example, we consider the algebraic 
spin liquid phase with staggered flux on a square lattice. 
We comment on the differences for the 7rFO and OFo 
phases. We do not consider the 7rFn phase, to avoid the 
complications related to the presence of an SU(2) gauge 
field. 



A. Framework 

We start with establishing the framework and intro- 
ducing the basic elements required to calculate the ultra- 
sound attenuation. These are the gauge field propagator 



FIG. 2. Double wavy line shows the gauge field propagator in 
the RPA approximation. Thin wavy line is the bare Maxwell 
propagator. 



and the phonon self-energy. 

The gauge field, strongly coupled to spinon emerges 
from a microscopic constraint and fluctuations around 
mean field anzatz. On a microscopic level it originates 
from constraint and is not dynamical. Non-trivial dy- 
namics of the gauge field is generated due to the cou- 
pling to fermions. 34,35 First, the Maxwell term will be 
generated while integrating out high energy degrees of 
freedom. Another contribution, which is singular com- 
pared to the Maxwell term, comes from the Dirac band 
touching and can be written via vacuum polarization op- 
erator for massless Dirac fermions. Total action of gauge 
field then becomes 



S n = 



1 



d 3 q 



a,(q)[(D M )-l(q) - ^{q^-q), 



where we use covariant notations in Euclidean space, q — 
{iuj^vfQ)- ^p,v{q) is the polarization bubble of Dirac 
fermions (see Fig. 2), and (D M )~l is the inverse Maxwell 
propagator of the gauge field (we work in the Lorentz 

gauge, k ■ a = 0), 



{D M )-m=(8, v - q ^)^(<i). (4(i) 



where II M (cf ) corresponds to the inverse propagator with- 
out tensor structure: 



n M («f) = 



(47) 



The polarization bubble at zero temperature (projector 
tensor structure is again omitted) is given by: 36-38 



(48) 



where we introduced the integer number of flavors of 
four-component Dirac fermions, N, in our theory. The 
physical case corresponds to N — 2 coming from spin. 
Action (45) translates into the total propagator for the 
gauge field, Fig. 2, given by 



D(q) = 



1 



N 1^2 

q 



2 /(Ne 2 ) 



(49) 



In what follows we will need the polarization bubble at 
finite temperature, which may be written as 



n, 



(50) 
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(51) 



(52) 



are orthogonal to each other and their sum reproduces 
the original zero-temperature tensor structure, Eq. (46). 
Explicit expressions for H A ' B along with detailed calcula- 
tions are available in the literature. 38 We need the asymp- 
totic expression for U A in the limit T ^> vp\q\ ^S> to: 



n 



27VTlog2 ( 1 . lc 



(53) 



Sound attenuation, a s , will be calculated from the self- 
energy of phonons, il p h(a;,q), arising due to interactions 
with spinons. More precisely, a s is given by the imagi- 
nary part of the retarded self- energy, 16 



" S = Im[n*(u;,|g|)] ( 



(54) 



with frequency and momentum related by the dispersion 
relation of the acoustic phonons, lo = v s \q\, where v s is 
the sound velocity. 

Let us discuss the approximations to be used in the 
calculation of the sound attenuation. For simplicity, we 
consider the clean case, i.e. we assume that the mean 
free path of spinons, I, is much larger than the ultra- 
sound wavelength, ql ^> 1. Also, we assume that the 
sound velocity is much smaller than the Fermi velocity, 
vp 3> v s . Under this condition, we immediately find 
that if lo and q axe the phonon energy and the wave 
vector, vpq = (vp/v s )oj ^> cu. Finally, in contrary to 
the case of spin liquid with a Fermi surface, 16 non-zero 
temperature is required to get non-vanishing sound at- 
tenuation in a Dirac spin liquid. This is a consequence of 
the energy and momentum conservation in the scattering 
process. Acoustic phonon cannot excite a particle-hole 
pair of spinons since the maximum momentum change 
for such a pair with energy u, Afc = uj/vf, is much 
smaller than phonon momentum, q — u)/v s . Therefore, 
we assume that the system is at a finite temperature 
T » (v F /v s )uj » ui. 

As noted above, there is a gauge field strongly coupled 
to the spinons. In order to have a control over its effects 
we artificially introduced the number of flavors, N, being 
equal to two in the physical case. Since gauge field prop- 
agator, Eq. (49), is proportional to l/N, effects of gauge 
field are suppressed for large N. We will perform calcu- 
lations of sound attenuation to the leading order within 
l/N expansion, commenting on the higher order terms. 



B. Sound attenuation in the sFn phase 

Having the basic ingredient for the calculation of sound 
attenuation at our disposal, we consider a s for longitudi- 
nal phonons in the sFn phase. As shown in Section III B, 



(a) - 




(b) 



oo 



FIG. 3. Contribution of spinons to the longitudinal sound 
attenuation. The bare contribution from spinons is given by 
the diagram (a). Diagram (b) accounts for the screening due 
to fluctuations of the gauge field. Black dots represent spinon- 
phonon interaction vertex, specified in the main text. 



there is no allowed coupling at the leading order in fc. All 
possible couplings at the next order are given by Eq. (33). 
For simplicity, we consider only first term in Eq. (33) [see 
Eq. (32)]. Combining Eqs. (23) and (32), corresponding 
spinon-phonon interaction vertex reads: 

Mi 1 \q)^~g^qMi 1 \q), 



(55a) 

M { k 1) (q) = qlk x r x +q 2 y k y ry, (55b) 

where q = (q x , q y ) is the unit vector pointing along q. In 
what follows, coupling constants with tilde are defined as: 



1 



.9 



2pUJ q 



9- 



(56) 



To leading order in 1/iV, the polarization operator of 
phonons due to interaction with spinons is given by the 
sum of two diagrams in Fig. 3 with the spinon-phonon in- 
teraction vertex from Eq. (55). Indeed, the first diagram 
in Fig. 3 has one fermionic bubble and is proportional to 
N . An extra fermionic bubble in the diagram Fig. 3 (b) is 
compensated by factor of 1 /N from the gauge field prop- 
agator. We have for the first contribution, Fig. 3 (a): 



a a = -$ 1} ]Vlmn*«( W ,9) 



(57) 



The imaginary part of the bubble diagram with spinons 
is calculated in the Appendix B and behaves as 
Imn fl W(aj, q) oc NujT 3 / (vpq) at the leading order. One 
can show that the contribution from the diagrams with 
an extra gauge field propagator, Fig. 3 (b), has the same 
order of magnitude as Fig. 3 (a). Thus, using disper- 
sion relation of acoustic phonons, we get the following 
estimate for the sound attenuation: 



u s Up 



(58) 



The value of the coupling constant, 9i , may be esti- 
mated from the sensitivity of the velocity of Dirac spinons 
to the changes of the lattice constant, a: 



9? 



dvp 
da 



v F . 



(59) 



Using this estimate, we obtain for the sound attenuation: 

2 



N 



T 



(60) 
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where defined as 



a 



(0) 



m ion v s VF 



(61) 



The Debye frequency has been estimated as u>n oc v s /a, 
and kx — T/vp is a wave vector of spinons with the 
energy equal to the temperature. The af \ introduced 
above, gives the estimate for the sound attenuation coef- 
ficient in the case of spinon Fermi surface if one substi- 
tute the Fermi momentum for the kx, kx — > kp. We see 
that in a Dirac spin liquid, contribution of spinons to the 
sound attenuation is suppressed compared to the Fermi 
surface case by two factors. The first factor, T//j,p <C 1 
is generic for any Dirac spin liquid and originates from 
the vanishing density of states at zero temperature in 
the Dirac spectrum. The second factor (T/ujpj) 2 , which 
is also expected to be smaller than one, arises due to 
peculiar form of spinon-phonon coupling. 

Finally, we comment on the next order in 1/N terms, 
contributing to the sound attenuation. There is a much 
larger number of diagrams at the order O(l). The most 
obvious are the vertex corrections, where gauge field 
dresses the interaction vertex of spinons with phonons or 
gauge field itself. Note, that if the spinon-phonon inter- 
action vertex corresponded to some conserved current, it 
would be protected from logarithmic corrections. 35 How- 
ever, this seems to be not the case here and in general 
we expect logarithmic corrections to arise at the order 
O(l). There is also another type of contribution O(l), 
which is more unusual. Indeed, in order to maintain the 
gauge invariance, k in Eq. (55) has to be extended to 
include the gauge field as well. This leads to the vertex 
where a phonon can generate a quanta of the gauge field 
in addition to the particle-hole pair of spinons. Similar 
type of vertex has been considered in Ref. 39. 



C. Sound attenuation in nF& and OFo phases 



attenuation, can be expressed via time component of elec- 
tron polarization bubble IIoo. In addition, one have to 
account for the effect of gauge field, including the scalar 
potential (unlike the case of spinons with Fermi surface, 
scalar potential of the gauge field is not screened by Dirac 
fermions). These two contributions to II p h are shown in 
Fig. 3, where now black dots correspond to the interac- 
tion vertex (62). Accounting for the both diagrams in 
Fig. 3, we get: 



n ph (g) 



[^ 0) ]V[noo(<?; 



no^iVOf) 11 ^?)], 

(63) 

where propagator and self-energy are taken on a phonon 
mass shell, uj = v s \q\. Using the finite-temperature ex- 
pression for the polarization operator, Eq. (50), we find 
that only H A contributes in the present case. Two terms 
in the sum in Eq. (63) partially cancel each other and we 
arrive at: 



n ph (<f ) = glq 



2 q* U M (q)U A (q) 

nM(q)_n^(g; 



(64) 



Since U A (q ) is more important than the Maxwell term, 
at the leading order we can neglect the latter term in 
the denominator, and get II p h (cf) oc — U M (q). Note, that 
this term is of order of 0(1), compared to the naive ex- 
pectation n p h((f) ~ O(N). Moreover, this term does 
not contribute to the imaginary part of the self-energy: 
Maxwell propagator originates from high-energy modes, 
whereas we are interested in the decay of phonons into 
low-energy Dirac-like spinons. Omitting the leading or- 
der term, and including next order contribution, we get 
the result 



n ph (g) = -glq 



2 q 2 [n M (<?) 



u A (q) 



(65) 



As we have shown above, the peculiar form of the cou- 
pling between phonons and spinons in the sFn phase 
leads to the suppression of the sound attenuation coeffi- 
cient by additional small factors. One may expect, that 
since in the 7rFx> and OFo phases longitudinal phonons 
couple to the density of spinons, the sound attenuation 
will be parametrically larger than for the sFn phase. Be- 
low, we are going to demonstrate that these naive expec- 
tations do not hold. Due to the presence of the gauge 
field, the density coupling gets screened and does not con- 
tribute to the sound attenuation at leading order in 1/N. 

Using explicit form of the coupling, Eq. (39), we write 
corresponding spinon-phonon interaction as: 

Mi \q)^g^qt. (62) 

The identity matrix corresponds to the spinon density, 
thus justifying the use of the term "density couping". 
The self-energy of phonons, required to find the sound 



which is proportional to 1/A^. Qualitatively, cancellation 
of two leading terms can be understood as an effect of 
screening due to gauge field. 

Now that we have shown that the contribution of the 
density coupling to the sound attenuation is proportional 
to 1/iV and thus negligible, we consider other terms in the 
coupling Hamiltonian, contributing at the order O(N). 
For the ttF& phase, these terms, listed in Eq. (42), are 
first order in k. Thus, the sound attenuation is expected 
to be of the same order as the results for the sFn phase, 
listed in Eq. (60). 

However, for the OFo phase there are couplings allowed 
without an extra k, see Eq. (44). Contribution from these 

couplings is expected to be of order of ai°^ [see Eq. (61)]. 
We note that contribution of gauge field vanishes in the 
present case. Indeed, it couples with an opposite sign in 
different valleys [note the presence of the extra [i z matrix 
in Eq. (44)], thus diagram in Fig. 3 (b) is identically zero. 



12 



V. DISCUSSION AND OUTLOOK 

We presented the general procedure for the derivation 
of the coupling between spinons and acoustic phonons in 
the Dirac spin liquid. Our procedure is based on the sym- 
metry arguments. Although general fermionic bilinear 
transforms under projective representation of the lattice 
symmetry group, spin singlet bilinears realize conven- 
tional (i.e., not projective) representation of the micro- 
scopic symmetry group. We found the decomposition of 
this representations into irreducible for 7r-flux and stag- 
gered flux phases on a square lattice, as well as for 7r-flux 
phase on kagome lattice and a Dirac spin liquid phase 
on a honeycomb lattice. By pairing corresponding irre- 
ducible representations with those for acoustic phonons, 
we were able to identify all symmetry allowed couplings. 
Note, that such decomposition can have other applica- 
tions. For instance, it can be used to derive symmetry 
allowed couplings to optical phonons or some other exci- 
tations. 

In a continuum limit all considered spin liquid phases 
have similar low energy Dirac excitations, and hardly 
can be distinguished. Nevertheless, the allowed inter- 
actions with phonons have different form. For the Dirac 
spin liquid phase on a honeycomb lattice the coupling 
to acoustic phonons is similar to the case of graphene. 
The only difference is that the coupling to the density of 
spinons, naively expected to be the largest, is screened 
by the gauge field (this is true for all Dirac spin liquid 
phases). As a result, for spin liquid phases on a square 
and kagome lattices considered in this work, the leading 
couplings contains an extra small parameter (T/ujd) 2 , 
compared to U(l) Dirac spin liquid on honeycomb lat- 
tice. Qualitatively, in these phases, the lattice deforma- 
tions with small wave vectors couple to the changes or 
anisotropics in Fermi velocity. Whereas in the case of 
zero-flux phase on honeycomb lattice such lattice defor- 
mations shift the position of Dirac points, acting similarly 
to the gauge field. 

The difference between couplings arises naturally from 
the fact that they are controlled by the representation 
of the corresponding symmetry group, acting on a lat- 
tice level. Thus the interaction of spinons with phonons 
retains some information about microscopic structure of 
the phase. It would be instructive to check if one can 
distinguish between different projective realizations of 
the same symmetry group by looking at couplings to 
spinons. The simplest example 40 of such two phases are 
two 1i 2 spin liquid phases on a square lattice (Z2A0013 
and Z2Azzl3 in notations of Refs. 23 and 41). This, how- 
ever, requires generalization of the present approach to 
the case of Z2 spin liquid phases, which is an interesting 
open question. Another open question is to understand 
the effect of projective realization of spin SU(2) symme- 
try, which has been proposed recently. 42 

In order to understand the perspectives of spinon- 
phonon interaction as a probe of fermionic spinons, we 
carried out a simple calculations within 1/N expansion. 



Assuming that our results can be extrapolated to the 
physical case N = 2, we see that the in a generic Dirac 
spin liquid exemplified by the zero flux phase in the 
honeycomb lattice, sound attenuation is suppressed due 
to vanishing density of states at zero temperature and 
as cx qT/ (mi on VFV s ). On the other hand, the pecu- 
liar form of spinon-phonon coupling in the 7r-flux and 
staggered flux phases contributes an additional suppres- 
sion of the form (T/lud) 2 . Nevertheless, the effect from 
phonons is still potentially observable, as the sound at- 
tenuation due to phonon-phonon scattering (caused by 
non-liearities) behaves as a ~ T 4 for T -C ojd- 

We note that our calculations should be viewed as 
a simple estimate due to the nature of approximations 
used. Currently, to the best of our knowledge there is 
no experimental data available on the sound attenuation 
in Dirac spin liquids. Provided such data becomes avail- 
able, more extensive theoretical work is required, in order 
to construct a realistic description. In particular, for the 
prospective spin liquid phase on a kagome lattice, 4,44 the 
clean limit, assuming mean free path I 3> q^ 1 does not 
apply. Also, the estimate for vp suggests that vp ~ v s , 
rather than vp 3> v s as was assumed. Another question, 
which can be relevant for a spin liquid on a kagome lattice 
is the effect of transition from U(l) to 1i 2 spin liquid and 
its possible manifestation in the ultrasound attenuation. 
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Appendix A: Elements of representation theory for 
relevant groups 

This Appendix provides background on the represen- 
tation theory, and gives more details for the symmetry 
groups used in the main text. It starts with a summary 
of the basic facts from the representation theory of fi- 
nite groups, which are extensively used throughout the 
paper. The reader interested in more details or deriva- 
tions of particular statements is referred to Refs. 45 and 
46. Next, the symmetry group of square and its exten- 
sion, relevant for the 7rFn and sFn phases, is consid- 
ered. Finally, the basic facts about the symmetry group 
of hexagon and the symmetry group of the 7rFO phase 
are discussed. 



1. Basic facts from representation theory 

We consider a point group Q, which contains hg ele- 
ments. Notion of conjugacy classes will be of great impor- 
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tance for us in what follows. Conjugacy class is defined 
as a complete set of mutually conjugate group elements, 
where two group elements g\ and g^ are defined to be 
conjugate if there exists another group element g%, such 
that g\ = g^ 1 o g 2 o <? 3 . In other words, if g belongs to a 
given conjugacy class, Ci, then for any group element 

Vffj-e^, gj 1 ogogj eCt, (Al) 

still is an clement from the conjugacy class Ci. Let us as- 
sume, that the group Q has ng conjugacy classes, denoted 
as C\,C2, ■■ .C ng . Each class contains N k elements, and, 
since each group element belongs to only one conjugacy 
class, we have J^fcii Nk — hg. The identity element, 
which is necessary present in any group is a conjugacy 
class itself, C\ = E = {1} and Nx = 1. For an abelian 
group, any element belongs to a separate conjugacy class, 
so that ng = hg, and Nx ,.... ng = L 

In what follows, our main interest will be in classi- 
fying representations of a given group. Representation 
of the group can be thought of as a mappings from 
the group elements to operators acting on some linear 
space, g — > R g which respects the group multiplication, 
R gi ■ Rg 2 = Rg l0g2 . If operators from a given repre- 
sentations cannot be represented as a direct sum of two 
operators acting on a smaller subspaces, this represen- 
tation is called irreducible. According to this definition, 
any representation D can be expressed as a direct sum 
of irreducible representations, 

D = ajflW © a 2 L> (2) 9 ... 9 a ng D^\ (A2) 

where non-negative integers a$ describe how many times 
a given irreducible representation is encountered in the 
decomposition. If D^> is not contained within D, corre- 
sponding aij is zero, a, = 0. In this way the problem of 
classifying all representations of a given group is reduced 
to a classification of all irreducible representations. 

The number of different irreducible representations for 
the group Q coincides with the number of its conjugacy 
classes, ng. Each irreducible representation, D^' is spec- 
ified by the value of its character for different conjugacy 
classes, defined as 

X {i) {C k ) =txR gCk , where g Ck eC k . (A3) 

According to the definition of the conjugacy class (Al), 
the value of X (^k) does not depend on the choice of a 
particular element gc k from the C k . Operators which act 
on a linear space can be expressed as matrices, and trace 
in Eq. (A3) is understood in this sense. 

Value of character for the conjugacy class which con- 
sists identity E = {1} is special, since it gives us the 
dimension of the corresponding irreducible representa- 
tion. To classify all irreducible representations of a given 
group, it is good to know not only the number of dif- 
ferent irreducible representations, but their dimensions 
as well. In such situation the following relation between 
the number of elements in the group, hg and the dimen- 
sions of all irreducible representations contained within 
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TABLE V. Irreducible representations of C4„ and their char- 
acters. 

the group, ti — x^ (E) turns out to be particularly use- 
ful: 

N c 

hg = Y J Ll (A4) 

i=l 

Typically only a few sets of integers {£x, ■ ■ ■ ,^jv g } satisfy 
this relation, and one can usually identify the correct set 
of dimensions by involving other considerations. 

Character table is a compact way of describing all ir- 
reducible representations of a given group. It is a square 
ng x ng table, where columns correspond to different 
conjugacy classes, and rows are labeled by different ir- 
reducible representations. The entry at an intersection 
of i-th row and j-th column is given by the value of the 
character for the i-th representation of the group ele- 
ments from the j-th conjugacy class. 

Using the characters table of a given group, one can 
easily find multiplicities Oj in the decomposition of a rep- 
resentation D into irreducible representations, Eq. (A2). 
Provided, characters of the representation D, x(Cfe), are 
known, we can find aj as 

, n e 

a^—J2N kX (l) *(C k )x(C k ), (A5) 
n ° k=i 

where hg is the number of elements in Q, and Nk is the 
number of elements in the corresponding conjugacy class. 

If representation D is obtained as the tensor product 
of two representations, let us say, E and F, D — E x F, 
the characters of D, x(Cfe) = x ExF (Ck), can be obtained 
as a product of characters for representations E and F, 

X ExF (C k )= X E (C k )x F (C k ). (A6) 

After this, one can easily apply formula (A5) to find the 
decomposition of the E x F into irreducible representa- 
tions. 



2. Group of square lattice and its representations 

Here we illustrate how the facts summarized above may 
be used to classify representations of the point symmetry 
group of square, Cx v and its extension, C' iv . 
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TABLE VI. Labeling of conjugacy classes of group C' 4v . Below each label, number of group elements, Nc, belonging to a given 
conjugacy class, as well as explicit form of these elements in Seitz notations are given. Vector a 3 is a short-hand notation for 
the sum of lattice vectors, a 3 = ai + 02. 



Point group of square Ca v and its representations 

We start with reviewing properties and representations 
of the point symmetry group of square, C^ v . This is the 
group of all symmetry operations, which leave square in- 
variant. It can be generated by rotations for tt/2 around 
the center of the square, and a reflection of x-axis, 

1Z X . 45 In total the group C 4v has hc 4v — 8 elements. In 
addition to rotations for angles multiple of tt/2, these in- 
clude reflections around x and y-axes, as well as !Z u ,v, 
standing for reflections relative to the planes containing 
vectors x ± y, TZ UlV = Tl Xt yll 7r /2- 

These elements can be split into total of nc 4ij = 5 
conjugacy classes. There are two conjugacy classes con- 
sisting of only one group element: trivial E = {1}, and 
C2 consisting of rotation for tt, C2 = {TZtt}- Each of the 
remaining three classes consists of two elements: C 4 = 

{^-7r/2j^37r/2}i C xy = V^xi Tty} , an d C uv — {7Z U ,1Z V }. 

Correspondingly, group C± v has five irreducible repre- 
sentations. Using Eq. (A4) we find that four of irre- 
ducible representations are one-dimensional and one is a 
two-dimensional. Characters of these irreducible repre- 
sentations are listed in Table V. One-dimensional repre- 
sentations are fully specified by their list of characters. 
Whereas two dimensional representation E\ corresponds 
to a transformation of a vector. If we denote the basis of 
E\ as (x,y), action of the group generators becomes 

TZ n /2 ■ x —> y, y — > — x, (A7a) 
1Z X : x —> —x, y —> y. (A7b) 

Using Table V we can easily find decomposition of the 
Kronecker product of E\ x E\ into irreducible represen- 
tations. 45 Only two non-zero characters of E\ x E\ are 
X ElxEl (E) = x ElxEl (C2) = 4. Now, using Eq. (A5) we 
can find that first four representations in Table V are 
contained once within E\ x E%: corresponding multiplic- 
ities are all equal to one, a, = 1/8 • (4 + 4) = 1. Whereas 
for Ei, we find corresponding a to be zero. This may be 
summarized as 

E x x E x = Ai®A 2 ®Bi®B 2 . (A8) 

Although formula (A5) gives us information about rep- 
resentations contained within E\ x E\, it does not give 



explicit expression for basis of these irreducible represen- 
tations. In the present case the explicit form of the basis 
may be easily guessed from physical arguments. Basis of 
each Ei in the product can be written as a two compo- 
nents of a vector, with the action of generators specified 
in Eq. (A7). Having components of two vectors (q x ,q y ) 
and (u x ,u y ), one can easily guess that the quantity, in- 
variant under all symmetries is the scalar product. Thus, 
q u = q x u x + q y u y is a basis of Ai component, contained 
in Eq. (A8). Basis for A 2 is also easy to guess, as it has 
to change sign under any reflections. Thus, it is given by 
the vector product, q x u = q x u y — q v u x . Finally, one 
can check that remaining combinations q x u x — q y u y and 
q x u y + q y u x realize the basis for Bi and B2 irreducible 
representations. This leads us to the Eq. (27) in the main 
text, which summarizes the above results. 

Group C 4v and its representations 

From the group C 4v we move to the group C' 4v — 
Gu/Git, which is the factor group of the space group 
of square lattice Grj over the group of translations for 
two unit cell vectors G 2 t ■ In other words, group C' 4v is 
defined as group C 4v with added translation operations 
t ai and t a2 . In order to specify this group we will use 
Seitz operators {R\t} defined as 

{R\t}r = Rr + t. (A9) 

The group 

G t = {{1|0}, {l|oi}, {t\a 2 }, {t\ai + a 2 }} (A10) 

is the subgroup of C' 4v and it contains h,Q t = 4 elements. 
The group C' 4v has hc t = h Gt -h 4v — 32 elements. It has 
no' — 14 conjugacy classes, which are listed in Table VI. 

Representations of C' 4v can be worked out us- 
ing the fact that it has a subgroup Gt- Conse- 
quently we can easily obtain five irreducible represen- 
tations, one-dimensional Ai^ and B12 along with two- 
dimensional Ei , from corresponding irreducible represen- 
tations of Ci v . For this we simply assume the action 
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TABLE VII. Irreducible representations of C' iv and their characters. The first eight representations are one-dimensional, the 
remaining six representations are two-dimensional. 



of translations to be trivial. Assuming that translation 
result in multiplying basis elements by minus one, we 
find additional four one-dimensional irreducible repre- 
sentations, denoted as A[ 2 and B[ 2 to emphasize that 
these are an extension of corresponding representations 
from Ci v . Analogous extension of E\ is denoted as E[. 
Remaining four two-dimensional irreducible representa- 
tions can be found explicitly using SU(A) generators 
given by {//, r\ //r- 7 } as a basis. Action of translations 
for representations E 2 ■ ■ ■ E 5 can be written as 

T x , y ■ x -» =fx, ±y, (All) 

so that T x T y — — 1. However, the transformation of basis 
under rotation and reflection are realized differently for 
each of these representations. For representations E 2 and 
E3 we have 

TZ n/2 : x^y, y^x, (A12a) 
1l x : x =px, y -> TV, (A12b) 

with the minus (plus) sign corresponding to E 2 (E 3 ). For 
Ei (E 5 ) we get: 

■ x -> TV, V -> ±x, (A13a) 
U x : x^tx, y^f ±y. (A13b) 

The character table may be easily calculated from here, 
and it is summarized in Table VII. 

From characters we determine the decomposition of 
different representations of C' iv on fermion bilinears into 
irreducible representations. Indeed, basis in the space 
of all possible fermion bilinears that are singlets in 
spin sector can be constructed using SU(4) generators 
, t 1 , [i l T*i } . Therefore this problem is equivalent to 
reducing adjoint representation of C' 4v on sixteen 4x4 
matrices {1, //, r l , /xV J '}. Representation is fully spec- 
ified by the action of generators. For the cases of the 



7rFn phase these are given by Eqs. (13)-(19). Whereas 
for the sFn phase, the reader is referred to Ref. 26. (Note, 
that there 1Z X is defined as a reflection with respect to 
the edge of square, whereas in our conventions reflection 
plane goes through the center of plaquette. Therefore 
1Z X from Ref. 26 coincides with T X 1Z X in our notations.) 
Calculating characters and applying Eq. (A5), we find 

=A 1 +A 2 + B[+B' 2 +E 1 +E[ 

+ E 2 +E 3 +E 4 + E 5 , (A14) 

for the 7rFn phase, where basis in terms of products of 
Pauli matrices for each irreducible componentis listed in 
Table I in the main text. We also obtained the same 
expressions for bases of different representations using 
the notations from Ref. 25. Analogously, for the 7rFn 
phase we have: 

G>t° - A 2 + B 1 + B\ + A' 2 + 2E[ + 2E 3 + E 4 + E 5 , 

(A15) 

with details on the basis listed in Table II. From here 
we immediately recover result of Refs. 25 and 26 that no 
invariant fermion bilinear terms exist in 7rFn and sFn 
phases. Indeed, G sF does not contain trivial represen- 
tation ^4i. Whereas, even though G vF contains At, as 
one can see from Table I it is not invariant under time 
reversal, T, nor under charge conjugation, C. 

3. Group of honeycomb and kagome lattices 

Since an extensive details for kagome and honeycomb 
lattices are available in the literature, 22,27 we only briefly 
summarize the basic facts for the symmetry group of the 
hexagon Cq v and its extension for the 7rF^ phase. More 
details for the honeycomb lattice can be found in Ref. 22. 
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Point group of hexagon 

For kagome and honeycomb lattices the relevant point 
group is that of a hexagon, denoted as Cq v . It has 
hc Bv — 12 elements and can be generated by the rota- 
tion TZ n / 3 and the reflection of y-axis, lZ y . It has six 
different conjugacy classes and six irreducible represen- 
tations, of which four are one-dimensional, and remaining 
are two-dimensional. Using characters of C 6v shown in 
Table VIII, we can write product of E\ x E\ as 

E 1 xE 1 = A 1 ®A 2 ®E 2 . (A16) 

Cg„ for kagome lattice 

Anzats for the algebraic spin liquid on Kagome lat- 
tice has a larger unit cell than the case without any 
fluxes. Thus, to classify fermionic bilinears, we again 
have to consider enlarged group, C' 6v , which is the G 6l) 



with added translations for primitive lattice vectors ai 
and ci2- 

The group C' ev (or, G S 2 in notations of Ref. 27) has 
been studied extensively and its conjugacy classes along 
with characters are listed in Tables III and IV in Ref. 27. 
Using this information, we may find the decomposition of 
the representation on bilinears as in Eq. (38) with bases of 
corresponding irreducible components listed in Table III. 

In the next order, we have to decompose the E\ x 
into irreducible representations. This leads us to 
Eq. (40) in the main text, where components A±, A2 
and E2, which are of interest for us originate from the 
tensor product of E\ with another E\, contained within 
Eq. (38) . This readily allows us to find the basis for these 
representations. 



Appendix B: Calculation of the polarization 
operator 

In this appendix we calculate the imaginary part of 
the polarization bubble. We work using assumptions, 
specified in the main text. In particular we restrict our- 
selves to the clean limit ql ^> 1 and assume the tem- 
perature to the the largest energy scale in the problem, 
T vpq ^S> cj- Note, that we use explicit value oi N = 2 
corresponding to spin. Since the polarization operator is 
proportional to N, one can easily restore the answer for 
the general case. 

We write the polarization operator, U^', correspond- 
ing to the interaction vertex M^ l '(q) as 



ImnW(iw n ,g) = 2TIm f (dk) ^ tt[M® (q)G k+q (iu m + ic^)M® q (q)G h (iu m )], (Bl) 

•* m 

where (dk) = dk x dk y / (2ir) 2 is the short-hand notation for the momentum integration measure. The interaction vertex 
M^'(q) as well as the Greens function are matrices in spinor space, and tracing in (Bl) goes over matrix indices. 
After analytical continuation, the imaginary part of the the matsubara sum of two Greens functions is written as, 

Im^[G fc+q (iw m + iuj n )] al3 [G k (iuj m )} 1 s = — — / dz Uanh ^= - tanh J Im[G£ +q (z + u)] a p Im[G^(z)] jS - 

(B2) 

where we restored internal indices. G k ' (z) stands for retarded (advanced) Greens function for real frequencies, 



Gl' A {z,k)= . *+Z r '* . (B3) 
(z ± 1O) 2 - Vpk 2 



In what follows, we will need the expression for the trace of numerators of two Greens functions with corresponding 
interaction vertices in Eq. (Bl). For the case of density coupling, defined in Eq (62), we have Mj^\q) = 1, and the 
trace is evaluated as: 

T (0) (z,w, k,q) = tr[l • (z + uj + v F r ■ (k + q)) ■ t ■ (z + v f t ■ fe)] = 4[(z + lu)z + (v F k) 2 + v 2 F kq cos 9], (B4) 

with 8 being the angle between vectors k and q. Note that there is an additional factor of two in (B4) from accounting 
for the (trivial) valley structure, whereas the factor of two originating from spin degrees of freedom is included in (Bl). 
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For the case of spinon-phonon coupling, arising in the next order of expansion in k, the m£ (q) is given by Eq. (55b) 
and the trace results in a cumbersome expression for TW(z,w, k,q), which will be not listed here. Using expression 
for the imaginary part of Green's functions, we have: 



Imn fl »(a;,g)=7rlm J (dk) J dz (tanh ~ - tanh T«(«, u, k, Q )-^ fe 



S(z + Vpk)8(uj + Upfc' — Upfe) 



S(z — v F k)S(uj — Vpk' + vpk) 5(z + Vpk)S(u> — vpk' — vpk) 5(z — v F k)S(uj + vpk' + vpk) 



vpk + uj 



vpk — UJ 



+ 



vpk + uj 



(B5) 



We drop last two terms in the square brackets since they correspond to interband transitions, and for oj <C Vpq they 
are not important. Also, we expand the difference between hyperbolic tangents, thus getting the derivative of the 
Fermi distribution function, denoted as n' F {z): 



lmU R(l) (uj,q) = 2ttuj (dk) / dzn' F (z)T^(z,u,k,q 



1 



Avpk 



S(z + v F k)S(uj + Vpk' — vpk) 
vpk — to 

S(z — vpk)S(uj — vpk' + Vpk) 
vpk + u> 



(B6) 



Using 5-functions, we may get rid of the integration over z. Integral over angle between vectors k and q, denoted as 
9, can be done using the following expression: 



vpk ± uj 



[ d96(±uj - v F \k + q\+ v F k)F(9) = 28(2k - g)-^- 
J vfikq\sm8$\ 



F(6t), 



where cos t 



2v F kq 



± 



v F q 



2k' 



(B7) 



This is valid in the limit when vpq 3> u). Note, that we included an extra factor 2 to account for two possible values 
of 6q (and 6q ), assuming that the F(9q ) is the same for both solutions. The integration over in Eq. (B6) yields: 



Imn fl 'W(w,q) 



Airvfiq 



dk 



3/2 



n' F (— Vpk) 



T^(-v F k,u,k,q)\ e=e 
k\ sin 8q\ 



-n' F (vFk) 



TM(v F k,w,k,q)\ g=6 
k\ sin 6+\ 



(B8) 



We notice, that expression in the square brackets in Eq. (B8) does not vanish if we put u) to zero within it for the 
case of density coupling [when T(z, k, 0) is given by Eq. (B4)]. In this case, accounting for the fact that n' F (v F k) for 
the vanishing chemical potential is an even function, we have: 



Imn fl (°»( W ,g) 



-iuj 



nv F q J q/2 



dkn' F {v F k)y/k 2 - (q/2) 2 



irvpq 



Tlog2. 



(B9) 



When calculating the integral we used the fact that the main contribution to the integral comes from v F k ~ T, thus 
we may neglect by q in the square root. This answer reproduces the results, available in the literature. 36_38 ' 47,48 
Recalling that this polarization operator is proportional to N, which was assumed to be N — 2 for this calculation, 
we reproduce the imaginary part of the result listed in the main text, Eq. (53). 

The calculation for the case of the next order coupling, M9~\q), requires more care. The answer depends on 
the direction of the phonon momentum, q. We define the <f> to be an angle of q relative to the £-axis, so that 
q = (cos 0, sin 0). Lengthy, but straightforward calculation gives for the polarization operator in this case 



Imn fl «(u;,g) 



uj sin 2(j> 
nq 



dk n'p(vpk) — 

q/2 FV V^ 2 ~ (9/2) 2 



2tt v 3 F q 



(BIO) 



Noteworthy, the answer is invariant under rotations of tt/2, as one may expect for our case. The angular dependence 
of (BIO) is very anisotropic, in particular, when q points along x or y-axes, the result vanish, indicating that the 
answer will be of higher order in uj. 
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